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Abstract. Functional determinants of differential operators play a prominent role 
in theoretical and mathematical physics, and in particular in quantum field theory. 
They are, however, difficult to compute in non-trivial cases. For one dimensional 
problems, a classical result of Gel'fand and Yaglom dramatically simplifies the 
problem so that the functional determinant can be computed without computing the 
spectrum of eigenvalues. Here I report recent progress in extending this approach to 
higher dimensions (i.e., functional determinants of partial differential operators), with 
applications in quantum field theory. 



1. Introduction : why are determinants important in QFT? 

This talk considers the question: what is the determinant of a partial differential 
operator, and how might one compute it? 

The physical motivation for this question comes from several areas of physics: 

• Effective actions, grand canonical potentials: These appear in relativistic and non- 
relativistic many-body theory and quantum field theory (QFT), and are generically 
given by an expression of the form Q = tr In G = In det G, where G is a Green's 
function, and the trace is over an infinite set of quantum states. 

• Tunneling and semiclassical physics: tunneling and nucleation processes can be 
analyzed in a semiclassical approach, wherein the leading exponential contribution 
involves the action evaluated on a classical solution, while the next-to-leading 
prefactor is the determinant of the fluctuation operator describing quantum 
fluctuations about the classical solution: 



• Gap equations: in many-body theory and relativistic QFT, one can find ground 
states by solving self-consistent Hartree-Fock or Schwinger-Dyson equations, 
varying the grand canonical potential or effective action with respect to a 
condensate field. Symbolically, the problem is of the form 
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where T> is some differential operator. Clearly, to be able to vary the determinant 
efficiently, one needs a clever way to evaluate the determinant. 

• Lattice gauge theory: integrating out quarks leads to fermion determinant factors 



which are needed for dynamical fermion computations in QCD. There are still 
fundamental challenges in evaluating such determinants within the Monte Carlo 
approach to lattice gauge theory, especially at finite density and temperature. 

• Faddeev-Popov determinants: the determinant of the Faddeev-Popov operator 
plays an important role in gauge fixing and confinement. 

• Mathematical physics: the determinant of a partial differential operator, such as a 
Dirac or Klein-Gordon operator, encodes interesting spectral information. 

2. What we know 

Exact results for functional determinants exist only for certain special cases. For 
example, the gauge theory effective action [the logarithm of the determinant of a 
Dirac or Klein-Gordon operator] for background fields with (covariantly) constant field 
strength can be expressed as simple proper-time integrals, building on the seminal work 
of Heisenberg and Euler [H [2], [3] . Mathematically, these expressions are closely related 
to the multiple zeta and gamma functions [I]. This constant curvature case is soluble 
because we can choose a Fock-Schwinger gauge with = — \F^ y x U) in which case 
the spectral problem decomposes into a set of harmonic oscillator problems, for which 
the spectral information is very simple. Similarly, in gravitational problems, for special 
spaces such as harmonic spaces or symmetric spaces, the spectral problem is sufficiently 
separable and simple that the determinant can be evaluated [5]. For example, the 
determinant of Dirac or Klein-Gordon operators on spheres, cones and tori are known 
in quite explicit form [5j [6] . Remarkably, this type of analysis can also be extended to 
general Riemann surfaces [3 [8] , with important applications to worldsheet problems in 
string theory and with beautiful applications in number theory and differential geometry. 

Things become more interesting when including general non-constant curvatures 
[either gauge or gravitational] . Then there are a few especially symmetric cases that are 
still soluble, but they are all one- dimensional. For example, in gauge theory, if the field 
strength points in a particular fixed direction, and has a magnitude that is a function of 
only one coordinate [say x^], and is of the form Fsech 2 (kx^), then the corresponding 
[Dirac or Klein-Gordon] spectral problem is hypergeometric, and the determinant can be 
computed explicitly. Also, for massless quarks in an instanton background, the fermion 
determinant is known. Such a computation was pioneered by 't Hooft for the single 
su(2) instanton [9], but has since been generalized to instantons of full ADHM form 
|10j . to finite temperature and to calorons with nontrivial holonomy [12] . These 
computations all rely on the fact that the Green's functions for massless particles in 
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an instanton background are known explicitly. The generalization to nonzero masses is 
nontrivial, and is discussed below in Section [HJ 

Some general bounds are known, based on lattice regularizations and generalizations 
of Kato's inequality for Dirac and Klein-Gordon operators. For example [T3l [b] 
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While very general, these bounds are still quite weak, and do not contain information 
about the phase of the determinant, which can be physically significant. Tighter bounds 
have been found recently in QED [15], but general results are still relatively limited. 

For Schrodinger operators with general potentials, or Dirac/Klein-Gordon operators 
with arbitrary gauge and/or gravitational backgrounds, we rely heavily on approximate 
methods. One such approximation is the heavy mass expansion. Define 

r °° d ^-e- m2s tr{e- sV } (5) 

Then an inverse mass expansion follows (after renormalization) from the small s 
asymptotic expansion of the heat kernel operator [16] 

1 
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Here the [D] are known functionals of the potentials appearing in T>. Similarly, a 
derivative expansion can be derived by expanding the heat kernel trace tr |e~ sI, | about 
the soluble constant background case, in powers of derivatives of the background. This 
corresponds to resumming all non- derivative terms in the inverse mass expansion. While 
very general, these expansions are asymptotic, with higher terms becoming rapidly 
unwieldy, and so have a somewhat limited range of application. 

Another approach is to define the functional determinant via a zeta function 



COO 
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A ^ 



(7) 



where the sum is over the spectrum of the relevant differential operator. Then, by the 
formal manipulations C'(0) = ~~ Sa In A = — In (J\x A), one defines the determinant as 



det 



-C'(o) 



(8) 



The problem becomes one of analytically continuing £(s) from the region in which it 
converges [typically dt(s) > d/2] to the neighbourhood of s = 0. For example, in a 
radially separable problem, one might estimate the spectrum using WKB phase shifts 
[and the corresponding spectral function p{k) = ^^f], thereby obtaining information 
about the zeta function. The zeta function approach is complementary to the so-called 
replica method, where one defines the logarithm of an operator by considering positive 
integer powers n of the operator, and then analytically continues to n = 0. There are 
also non-trivial interesting one dimensional quantum mechanical problems for which the 
determinant can be computed from the zeta function [T7j . 
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There are several special features of one- dimensional problems that deserve mention. 
First, the issue of renormalization does not really enter in Id. Second, the theory 
of inverse scattering allows us to characterize a potential in terms of the associated 
scattering data. This is especially useful for solving variational and gap equations, as 
the variational problem can then be reformulated in terms of the scattering data, which 
is more direct. Beautiful applications of this idea occur in interacting 1+1 dimensional 
QFT's such as the Gross-Neveu model, Sine-Gordon model and their generalizations 
[T8l [T9] . [Inverse scattering has been generalized in higher dimensions for radial problems 
[20J]. A third special feature of Id problems is the Gel'fand Yaglom theorem [2TJ [22] . 
which provides a way to compute the determinant of an ordinary differential operator 
without computing its eigenvalues. The basic idea is that in Id one can define the 
Green's function without using an eigenfunction expansion, simply from the product of 
two independent solutions. 



3. Gel'fand- Yaglom theorem: ordinary differential operators 



The simplest form of this result can be stated as follows [211 [22J [231 EH 125], [26] . Consider 
a Schrodinger operator on the interval x G [0, L], with Dirichlet boundary conditions: 

d 2 



dx 2 



V(x) 



ijj{x) = Xip(x) 



(9) 



Then to compute the determinant it is in fact not necessary to compute the infinite 
discrete set of eigenvalues {Ai, A2, • • •}• Instead, solve the related initial value problem 
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Then 
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It is worth pausing to appreciate the simplicity of this result, as well as its practical 
utility. It is straightforward to implement numerically. The determinant follows, without 
having to compute any eigenvalues, let alone multiply them all together! 

Actually, (fTTj) is not quite right, as we can only define the ratio of two determinants 
[27], and (TTTT) should be understood in this sense. Often in physical applications one 
is computing the determinant relative to the corresponding determinant for the free 
operator, for which the initial value problem is usually known in closed form [in this 
case it is numerically better to solve directly the initial value problem for the ratio] . As 
an illustration, consider the following simple example of the massive Helmholtz operator 
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The Dirichlet spectrum is A n = m 2 + (jjf\ , so from the eigenvalues: 
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On the other hand, if we use the Gel'fand-Yaglom result, we solve the initial value 



problems: 



for the massless operator. Then ( 1TTT) implies 



0, to find 4>(x) 



and -^0o 



0, to find 4>o(x) = x 
(13) 



These clearly agree, but the point is that if m 2 were replaced by a nontrivial potential 
V(x), the first approach, from the eigenvalues, would be extremely difficult, while the 
Gel'fand-Yaglom approach is still easy As another example, consider the Poschl- 
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Teller system with H = [- 



dx 2 



m 



j(j + l)sech 2 x], with integer j. This is a 



reflectionless potential with j bound states at E 
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shift S(k) = 2J2j=i arctan(//fc). The determinant follows directly from the spectrum as 
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On the other hand, the solutions to the Gel'fand-Yaglom initial value problems (fTUj) are 
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on the interval x G [—L,L]. The P m and are Legendre functions. Then 
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efaa i ... j - ~r(j=0)(L) V{m - j)T(m + j + 1) 

Similarly, one can define a general reflectionless potential with iV bound states at 



2 In det A(x) where is the iV x iV 



E = m 2 — k 2 , (i = 1, ... TV), by V(z) = m 2 
matrix A™ = <5j,- + ^ e -(«»+%-) a: . The c* are (constant) moduli which affect the shape 

of the potential, but not the spectrum. From the phase shift 5(k) = 2 J2iL\ arctan(Kj/fc), 
one finds the determinant 
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It is instructive to verify this numerically for various potentials using fllOp and (11 II) . 

The basic result (Hip generalizes in several ways. First, it generalizes to other 
boundary conditions [2H EB] , characterized by two 2x2 matrices M and N as 
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Then the infinite dimensional determinant reduces to a simple 2x2 determinant 
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For example, 

Dirichlet : M = 

Neumann : M 



V(x) 
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(20) 
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det 



M (2 )(L) 

u' (1) (L) 



(n (1) (L) + n' (2) (L)) 
(u ( i)(L) +u' (2) (L)) 



Periodic (P) : M = 1, N = -1 det = 2 

Antiperiodic (AP) : M = 1, N = 1 => det = 2 

In the P/AP cases we recognize («(!)(L) +tt'^(L)) as the Floquet-Bloch discriminant, 
which takes values ±2 for P/AP eigenfunctions. Second, these results generalize 
straightforwardly to coupled systems of ODE's [2H [28], and also to linear ODE's of 
any order (29]. Third, they generalize to Sturm-Liouville problems [28] • For example, 
consider the radial operator in d dim., with Dirichlet b.c.'s 



M 



d 2 



(0 



l + ^)(l + ^ 



d-l 



V{r) 



dr 2 r 2 

Define the function 0m (r) as the solution to the radial initial value problem 
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Then (this should also be understood as applying to ratios of determinants) 
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As an explicit example, consider the 2d radial Helmholtz problem, with eigenvalues in 



terms of the zeros j^ )n of the J; Bessel function: Am^. 
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On the other hand, if we use the Gel'fand-Yaglom result (|23|) . we readily solve the 
relevant initial value problems to find 4>(i) i n terms of the modified Bessel function If. 

free / r>\ / „\ I \ I 



The agreement between fl24]) and fl25|) expresses the product formula for Ii(x) 



4. Gel'fand-Yaglom theorem: partial differential operators? 

In generalizing from ordinary to partial differential operators, it is natural to consider 
first the case of separable operators. But the naive extension fails. Formally, we simply 
take a product over the radial determinants [each of which is evaluated as above] for all 
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possible angular momenta, with the appropriate degeneracies. However, this product 
diverges [23] . For example, for the 2d radial Helmholtz system ( 125]) 
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The degeneracy factor is 1 for / = 0, and 2 for all I > 1, so the angular momentum 
sum over / diverges. This should not be so surprising from a physical point of view, 
as we expect to require renormalization in dimension higher than 1. To formulate the 
problem more precisely, consider the radially separable operators J\A = —A + V(r), and 
J\A hcc = —A, where A is the Laplace operator in IR^, and V(r) is a radial potential 
vanishing at infinity as r~ 2 ~ e for d = 2 and d = 3, and as r~ 4 ~ e for d — 4. Since 
V = V(r), we can separate variables and consider the Schrodinger-like radial operator 
My) in (l2TI) . For dimension d > 2, the radial eigenfunctions ^(z) have degeneracy 

(2/ + d-2)(/ + d-3)! 
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Each term in the sum is computed straightforwardly using ff23]) . but the I sum is 
divergent. However, this divergence can be understood physically, leading to a finite 
and renormalized determinant ratio [301 : 
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Here 7 is Euler's constant, and /x is a renormalization scale, which is essential for physical 
applications. A conventional renormalization choice is to take \i = m in (129]) . In each 
of ( 129]) . the sum over I is convergent once the indicated subtractions are made. The 
function 4>(i)(r) is defined in ( 122]) . 

Notice that the results ( 129]) state once again that the determinant is determined by 
the boundary values of solutions of [M + m 2 ] = 0, with the only additional information 
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being a finite number of integrals involving the potential V(r). We also stress the 
computational simplicity of fl29|) ; no phase shifts or eigenvalues need be computed to 
evaluate the renormalized determinant. 

These results have been derived [30] using the zeta function formalism. Scattering 
theory [31] and standard contour manipulations [281 132] yield an integral representation 
[valid for 3?(s) > d/2] of the zeta function in terms of the Jost function fi(ik): 

OO 00 Q 

C(s ) = ^ deg(/; d) J dk [k 2 - m 2 ]" s ^ In /,(<*) . (30) 

The technical problem is the analytic continuation of (1301 to a neighborhood about 
s = 0. This analytic continuation relies on the uniform asymptotic behavior of the Jost 
function fi(ik), which follows from standard results in scattering theory [31] : 
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ln/,(iJfe) = / drrV{r)K v (kr)I v (kr) - J drrV{r)K 2 {kr) J dr' r' V(r')I 2 (kr') + 0(V 3 ) 



where v = I + | — 1. Then the calculation of the asymptotics of the Jost function 
reduces to the known uniform asymptotics of the modified Bessel functions K v and I v . 
Using these asymptotics, we can define In fi sym {ik) as the 0(V) and 0(V 2 ) parts of 
this uniform asymptotic expansion, and then by construction, 

oo 

C/(0)= - "£deg(l- d) [In Mim) - In fr m (im)} . (31) 

1=0 

is now convergent. Adding back the subtracted terms in a dimensional regularization 
scheme, leads EED to the results in 



5. Comparison With Feynman Diagram Approach 



These determinants can also be derived in a Feynman diagrammatic approach, where 
the perturbative expansion in powers of the potential V is 
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Here the solid dots denote insertions of the potential V. If the potential is radial, then 
with dimensional regularization one can express the renormalized 4d determinant as [33] 
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where the finite contributions from the first and second order Feynman diagrams in the 
MS scheme are [33] (note the small typo in equation (4.32) of [33]): 
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Here V"(g) is the 4d Fourier transform of the radial potential V(r). These expressions 
( 1331) and ( 1341) look quite different from ( 1291) . but in fact they can be shown to be 
completely equivalent, using some Bessel identities [30]. However, it is clear that (1291) is 
much easier to evaluate, especially if (as often happens) V(r) is known only numerically. 
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6. Removing Zero Modes 



In certain quantum field theory applications the determinant may have zero modes, 
and correspondingly one is actually interested in computing the determinant with 
these zero modes removed, together with a collective coordinate factor [23]. For 
example, consider the case of a self-interacting scalar field theory in d dim Euclidean 
space with Euclidean action S[<&] = J d d x [0 M <3>0 M $ + [/($)], where there is a [radially 
symmetric] classical solution $ c j(r) solving A$ = with $'(0) = 0, and $ — > as 
r — > oo. The fluctuation operator about this classical solution is a radial operator, with 
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In the / = 1 sector the corresponding radial operator 



Ai(i=i) has a d-fold degenerate zero mode associated with translational invariance: 

^W = ^oj ' (35) 

Including the collective coordinate contribution [23] due to translational invariance, the 
required determinant factor is [the prime denotes exclusion of zero mode(s)] 
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Since the determinant is expressed in terms of a solution to an ODE, we can use simple 
ODE theory to show that (generalizing the d — 1 [34] and d — 4 [35] analyses) 
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This expresses the zero mode factor solely in terms of the asymptotic properties of the 
classical solution $ c /, which are already known numerically as part of the determination 
of $ c z and the fluctuation potential V(r). No further computation is needed. 
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7. Application 1: fluctuation determinant for false vacuum decay 

The phenomenon of nucleation drives first order phase transitions in many applications 
in physics. The semiclassical analysis of the rate of such a nucleation process was 
pioneered by Langer [36] , who identified a semiclassical saddle point solution that gives 
the dominant exponential contribution to the rate, with a prefactor to the exponential 
given by the quantum fluctuations about this classical solution. The nucleation rate is 
given by the quantum mechanical rate of decay of a metastable "false" vacuum, into 
the "true" vacuum, $ + . Decay proceeds by the nucleation of expanding bubbles of true 
vacuum within the metastable false vacuum [36], [37J EH EHJ EUJ E]. The semiclassical 
prefactor requires computing the determinant of the differential operator associated 
with quantum fluctuations about the classical solution. This is a technically difficult 
problem, but it is ideally suited to the above results for determinants of radially separable 
operators, because the fluctuation operator is radial. The expression fj4*2]) below can be 
applied to any metastable field potential U(Q) [35], without relying on the "thin-wall" 
limit (the limit of degenerate vacua) . For definiteness, consider the specific (rescaled) 
quartic potential: U(&) = \<& 2 — |$ 3 + f$ 4 , where a characterizes the shape of the 
potential [the thin- wall limit is a — > 1]. The decay rate per unit volume and unit time 
is 



where the prime on the determinant means that the zero modes (corresponding to 
translational invariance) are removed. Here $ c i is a classical solution known as the 
"bounce" solution [39|, defined below, and the prefactor terms in ( |39|) correspond to 
quantum fluctuations about this bounce solution. The second term in the exponent, 
<5rfiS[$ c i], denotes the counterterms needed for renormalizing the classical action S c \. 
The first step in computing the false vacuum decay rate 7 is to find the classical bounce 
solution, § c \(r), which is a radially symmetric stationary point of the classical Euclidean 
action, interpolating between the false and true vacuum. The bounce $ c i(^) solves the 
nonlinear ordinary differential equation 



with boundary conditions = O5 &ci(r) — ► = 0, as r — > 00. <&d(r) must be 

computed numerically, for example by the method of shooting, adjusting the initial 
value $ c i(0) until the boundary conditions at r = 00 are satisfied. Given the bounce 
solution, & c \(r), the corresponding radial fluctuation potential is 



There are three different types of eigenvalue of the fluctuation operator, each having 
a different role physically and mathematically, (i) Negative Mode: The I = sector 
has a negative eigenvalue mode of the fluctuation operator, and is responsible for the 
instability leading to decay. This mode contributes a factor to the decay rate 7 related 
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(41) 
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to the absolute value of the determinant of the I = fluctuation operator 
(ii) Zero Modes: In the I = 1 sector, there is a four-fold degenerate zero eigenvalue 
of the fluctuation operator. Integrating over the corresponding collective coordinates 
produces the factors of |^ in ([39]) . In computing the rate 7, we need the determinant of 
the fluctuation operator with the zero mode removed [36 1 139 1 14*0] . This can be evaluated 
simply using the formula fl38l) . in terms of the asymptotic behavior of the classical bounce 
solution $ c i(r). (iii) Positive Modes: For I > 2, the fluctuation operator has positive 
eigenvalues, each of degeneracy (I + l) 2 . For each I, the associated radial determinant is 
computed numerically using (123]) . But the sum over I is divergent and the answer must 
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Figure 1. Comparison as a — > 1 of the exact numerical result 
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be renormalized. Using the MS scheme, the renormalized effective action is [35J: 
1 



r MS = g ln | T (0)(°°) 



2 In 
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ln^ 
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(42) 



--/ drrV(r) + — drr 3 V(V + 2)(- 
4 Jo 16 Jo \2 

The first term is from the negative mode, the second from the zero mode(s) [using (1381) ] . 

the third term from the numerical computation of the higher modes, and the last term 

is the renormalization counterterm contribution. The expression (I42p can be applied to 

any field potential U(<&), without relying on the thin-wall limit. A comparison with the 

analytic thin-wall computation of Konoplich and Rubin [12] is plotted in Fig. [TJ and 

shows excellent agreement. This approach can also be extended to false vacuum decay 

in curved space, where the structure of the bounce solutions [33] , and the corresponding 

fluctuation equations [44], is much richer. 
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8. Application 2: determinant for massive quarks in instanton background 

The fermion determinant in an instanton background is known in the massless (chiral) 
limit j9], and in the heavy quark limit [15]. Here we compute it for any m [16]. The 
relevant part of the renormalized effective action, rf en (m), behaves as 

«(l/2) + - (lnm + 7 - ln2) m 2 + . . . , m — > 

r ren(ra) = \ \ nm 1 17 2 32 7916 ^ 

H a ; o + ■ • • , m — > oo 



6 75m 2 735m 4 2835m 6 148225m 8 
where a(l/2) =~ 0.145873. The small mass expansion [9] is based on the known 
massless propagators in an instanton background, while the large mass expansion [35] 
follows most easily from the Schwinger-DeWitt or heat kernel expansion. 

Since the instanton is self-dual, we can simplify things slightly and work with 
the scalar Klein-Gordon operator. For a single su{2) instanton, the spectral problem 
separates into partial waves [9], with radial Schrodinger-like operators (take isospin |) 



H 



d 2 3 d 41(1 + 1) , A(j + l + 3 



(44) 



dr 2 r dr r 2 r 2 + 1 (r 2 + 1) 

Here I = 0, |, 1, |, ■ • • , and j = |/±||, and there is a degeneracy factor of (2/ + l)(2j + 1). 
Since the "potential" involves I and j, we cannot directly use the result (|2"9"1) . In order 
to compare with 't Hooft's massless quark analysis [9], we use a Pauli-Villars regulator 
A, and define the renomalized effective action by charge renormalization. To extract 
the renormalized effective action we need to consider the A — >• oo limit in conjunction 
with the infinite sum over /. Split the partial wave sum into two parts: 

L oo 

rt(m)= J2 r A >(0 H+ E r i(oH ( 45 ) 

1=0,1... l=L+± 

where L is a large but finite integer. The first sum involves low partial wave modes, 
and the second sum involves the high partial wave modes. For the low partial wave 
modes, we can remove the regulator and evaluate the (finite) sum using the (numerical) 
Gel'fand-Yaglom result (!23|) . But this sum diverges quadratically with L. In the second 
sum in ( l4"5l) we cannot take the large L and large A limits blindly, as each leads to a 
divergence. Radial WKB is a good approximation for the high / modes. With WKB we 
can compute analytically the large A and large L divergences of the second sum in ( 1451) . 
using the WKB approximation for the corresponding determinants 

OO 1 

r£ i(0 (A;m)~-lnA + 2L 2 + 4L 



l=L+i 



+ 




127 1, m 2 2 , m 2 , 

In 2 H m In 2 H In m 

72 3 2 2 



+ 0[j) (46) 



It is important to identify the physical role of the various terms in (146]) . The first term 
is the expected logarithmic charge renormalization counterterm. The next three terms 
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Figure 2. Plot of our numerical results for ff cn (m) from (|47p. compared with the 
analytic extreme small and large mass limits [dashed curves] from (T45)) . The dots 
denote numerical data points from (|47p . and the solid line is a fit through these points. 



give quadratic, linear and logarithmic divergences in L. These divergences cancel exactly 
against corresponding divergences in the first sum in (|4"5j) . which are found numerically. 

Combining the numerical results for low partial wave modes with the radial WKB 
results for the high partial wave modes the renormalized effective action rf en (m) is 




This is finite for any mass m. Figure [2] shows these results for ff en (m), compared with 
the analytic small and large mass expansions in (1431) . The agreement is remarkable 
|46j . It is a highly nontrivial check on the WKB computation that the correct mass 
dependence rises to interpolate between the large mass and small mass regimes. As an 
interesting analytic check, the formula fT4T|) also provides a simple computation [16] of 
't Hooft's leading small mass result a(l/2) = -f| - |ln2 + \ - 2C'(-1) = 0.145873.... 

9. Concluding remarks 

In this talk I have reviewed some of the reasons for considering determinants of partial 
differential operators in quantum field theory. Mathematically, much is known for 
ordinary differential operators, especially through the work of Gel'fand-Yaglom [2TJ [22J 
and its extensions [241 [26l [28] . but considerably less is known for partial differential 
operators with non-trivial potentials. The new mathematical results [30] reported here 
are the formulas in (T29]) . which provide simple expressions for the determinant of a 
radially separable partial differential operator of the form —A + m 2 + V(r). This 
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generalizes the (Dirichlet) Gel'fand-Yaglom result (II ip to higher dimensions. These 
results lead to many direct applications in quantum field theory, where they extend 
the class of solvable fluctuation determinant problems away from the restrictive class of 
constant background fields, or one dimensional background fields, to the more general 
class of separable higher dimensional background fields. This includes, for example, 
applications in quantum field theory to the study of quantum fluctuations in the presence 
of vortices, monopoles, sphalerons, instantons, domain walls (branes), etc ... 

A number of further generalizations could be made. First, it is practically important 
to include directly the matrix structure that arises from Dirac-like differential operators 
and from non-abelian gauge degrees of freedom. This has been done in various ways for 
certain applications [I7J HH HH EH [51], but a simpler unified formalism along the lines 
of the explicit formulas ( !29l) is probably possible, as in one dimension [52]. Second, 
a practical extension to finite temperature and density would be useful in various 
physical applications. Third, the biggest technical challenge is to ask if the restriction 
of separability can be relaxed. A promising current approach is the numerical worldline 
loop method [53]. Without symmetry, computing functional determinants analytically 
remains a difficult problem, and any progress will be interesting. 
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